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In order to clarify exotic multipole properties of filled skutterudites, we evaluate multipole suscep- 
tibility for n=l~5, where n is the local /-electron number, on the basis of a multiorbital Anderson 
model constructed using the j-j coupling scheme. For n=l, magnetic fluctuations dominate over low- 
temperature electronic properties, while for n=2 and 4, electronic states are dominated by both magnetic 
and quadrupole fluctuations. For n=3 and 5, octupole fluctuations are found to be significant, depending 
on the crystalline electric field potential. We discuss possible relevance of the results to actual materials. 

KEYWORDS: Multipole susceptibility, Filled skutterudites, Multiorbital Anderson model, j-j coupling 
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Recent discovery of heavy-fermion superconductivity in 
PrOs4Sbi2 has triggered a rapid increase of research activi- 
ties on filled skutterudite compounds. 1 ' 2 This material is be- 
lieved to belong to the group of strongly correlated electron 
systems, but the superconductivity cannot be simply under- 
stood in terms of d-wave pairing mediated by antiferromag- 
netic spin fluctuations, since the inverse of the nuclear spin- 
lattice relaxation time for PrOs4Sbi2 does not show the T 3 
behavior at low temperature T characteristic of a line-node 
gap. 3 In addition, PrOs4Sbi2 exhibits several exotic supercon- 
ducting features such as multiple superconducting phases, 4 
point-node behavior in the gap function, 4 5 and the breaking 
of time-reversal symmetry as detected by /iSR experiments. 6 

It is premature to settle on the mechanism of supercon- 
ductivity in PrOs4Sbi2, but recently, the possibility of exotic 
pairing due to quadrupole fluctuations has been suggested ex- 
perimentally 7 and theoretically. 8 In fact, neutron scattering 
experiments have revealed that a quadrupole ordered state 
is induced by a magnetic field in the vicinity of the super- 
conducting phase of PrOs4Sbi2, 9 implying a potential role 
for quadrupole fluctuations in superconducting pair forma- 
tion. Moreover, PrFe4Pi2 exhibits a second-order phase tran- 
sition at 6.5K, 10 but this transition is considered to be due to 
antiferro- quadrupolar ordering from neutron scattering ex- 
periments 11 and theoretical efforts. 1213 It has been gradually 
recognized that the quadrupolar degree of freedom plays a 
crucial role in filled skutterudite compounds. 

From the theoretical viewpoint, multipole phenomena in /- 
electron systems have been investigated phenomenologically 
on the basis of the LS coupling scheme. This Heisenberg- 
like model for relevant multipole moments could successfully 
explain some of the experimental results for filled skutteru- 
dites. In particular, the phase diagram for PrOs4Sbi2 has been 
reproduced 14 in a quasi-quartet system with Ti ground and 
excited states with a very small excitation energy 15 for a 
crystalline electric field (CEF) potential with 21 symmetry. 16 
However, the microscopic approach to electronic properties 
of filled skutterudites has been very limited mainly due to the 
complexity of multipole degrees of freedom. 

In order to overcome such difficulties, it has been pro- 
posed to construct a microscopic model for /-electron sys- 
tems based on the j-j coupling scheme. 17 With this approach, 
for instance, the microscopic origin of octupole ordering in 



NpC>2 has recently been established. 18 Furthermore, the mag- 
netism and superconductivity of filled skutterudite materials 
have been studied from a microscopic viewpoint on the ba- 
sis of a multiorbital Anderson model. 19 ' 20 However, the effect 
of multipole fluctuations has not been investigated satisfacto- 
rily at the microscopic level, even though their potential role 
has been emphasized in research on filled skutterudites. In the 
j-j coupling scheme, we can microscopically evaluate phys- 
ical quantities related to multipoles by treating them as com- 
bined spin and orbital degrees of freedom. Such a study is be- 
lieved to open a door which leads to a new stage for multipole 
physics of /-electron systems. 

In this Letter, we attempt to clarify what kind of multipole 
fluctuations are dominant in filled skutterudites, by analyz- 
ing a multiorbital Anderson model based on the j-j coupling 
scheme for n=l~5, where n denotes local /-electron num- 
ber. For this purpose, we evaluate susceptibilities for fifteen 
kinds of dipole, quadrupole, and octupole moments up to rank 
3 using a numerical technique. In this way, we find that for 
n=l, magnetic fluctuations are dominant at low temperatures, 
while for n=2 and 4, electronic states are dominated by mag- 
netic and quadrupole fluctuations at low temperatures. For 
n=3 and 5, octupole fluctuations are found to appear for some 
ranges of CEF parameters. We provide some explanations and 
predictions on filled skutterudites from the present results. 

In general, in the j-j coupling scheme, a large spin-orbit 
interaction is first assumed, and then, only the j=5/2 sextet 
is taken into account. Concerning the applicability of the j-j 
coupling scheme to rare-earth compounds, readers can con- 
sult Ref. 20. Note that the main conduction band of filled skut- 
terudites is a u with xyz symmetry, 21 where the subscript "u" 
means ungerade, i.e., odd parity for space inversion. Then, the 
Anderson Hamiltonian is written as 



H = 



ka 



EfeC^Cfeo- + y^CVc^/cq + h.c.) + H\ oc 



(1) 



ka 



where Sk is the dispersion of a u conduction electrons with Tj 
symmetry in terms of the j-j coupling scheme, is the an- 
nihilation operator for conduction electrons with momentum 
k and spin a, / 7cr is the annihilation operator of / electrons 
on the impurity site with pseudospin a and orbital 7, V is 
the hybridization between conduction and / electrons with T7 
symmetry, and H\ oc is the local /-electron interaction term. 



2 J. Phys. Soc. Jpn. 



Letter 



Takashi HOTTA 



The orbital index 7 is introduced to distinguish three kinds 
of the Kramers doublets, two T 8 and one T 7 . Here "a" and 
"b" denote the two r 8 's and "c" indicates the T 7 . Throughout 
this paper, we set V=0.05, where the energy unit is taken as 
D, half of the bandwidth of the conduction band. From band- 
structure calculations, the bandwidth is estimated as 2.7 eV in 
PrRu 4 Pi2, 22 indicating that D=1.35 eV. 

The local /-electron term H\ oc is given by 17 

ffioc = 5> 7 /t ff / 7CT + (i/2) £ E^v™, 

X f fieri f-y2<T2 f fjiCT! I (2) 

where _B 7 is the CEF potential, which is already diagonalized. 
It is convenient to introduce a level splitting A between r 7 
and r 8 as A=£>r 8 — B^ 7 - The Coulomb integral U in the j-j 
coupling scheme is expressed in terms of three Racah param- 
eters, Eo, Ei, and E^} 1 To set the local /-electron number, 
we adjust the /-electron chemical potential for each n. 

In order to reproduce the CEF scheme for Pr-based filled 
skutterudites, we prefer positive A. In particular, the quasi- 
quartet situation with singlet ground and triplet excited states 
is reproduced by taking a small positive A. However, since the 
difference between (\ and Th is not fully included in the j-j 

coupling scheme, 20 the excited state for n=2 is always T4 for 

(2) 

A>0, although T\ ' (a mixture of T 4 and T 5 ) is the ground 
state in actual material. The extent of the mixture explains 
the difference between PrOs4Sbi2 and PrFe4Pi2, 23 but here 
we simply ignore such effect. It is one of future problems to 
include the effect of Th in the j-j coupling scheme. 

In order to discuss multipole properties, we evaluate multi- 
pole susceptibilities for / electrons, given by 

where Ar 7 is the multipole operator, X denoting the multi- 
pole symbol, T is the irreducible representation with 7 to dis- 
tinguish degenerate representations, E n is the eigenenergy for 
the n-th eigenstate |n), and Z is the partition function given 
by Z=12n e~ B "/ T . For j=5/2, we can define multipole op- 
erators up to rank 5 in general, but we are primarily inter- 
ested in multipole properties from the T g quartet. Thus, we 
concentrate on multipole moments up to rank 3. 24 Note that 
in the j-j coupling scheme, physical quantities for multi-/- 
electron systems can be evaluated from one-particle operators 
on the basis of quantum-field theory, as we have performed 
for d-electron systems. We have the clear advantage that it is 
not necessary to redefine Ar 7 depending on the ground-state 
multiplet, when we change local /-electron number. 

Now let us show explicit forms for multipole operators fol- 
lowing Ref. 24. As for dipole moments with T 4u symmetry, 
we express the operators as 

J4ux Jxi J^uy — Jyi J^uz — Jzi (4) 

where J x , J y , and J z are three angular momentum opera- 
tors for j=5/2, respectively. Concerning quadrupole moments, 
they are classified into r 3g and r 5g , where g denotes gerade. 
We express the r 3g quadrupole operators as 

Oz g u = (2J 2 — J x — J y )/2, ,~ 
3gv = V3(JJ - J 2 )/2. W 



For the r 5g quadrupole, we have the three operators 

05g£ = \/3 JyJ Z / 2, 

05 g 7j = \/3 J z J x /2, (6) 

05gC = V3 JxJy/%, 

where the bar denotes the operation of taking all possible per- 
mutations in terms of cartesian components. 

Finally, regarding octupole moments, there are three types 
as r 2u , T4 U , and T 5u . We express the r 2u octupole as 

T 2u - Vl5J x J y J z /6. (7) 

For the T4 U octupole, we express the operators as 

T/luX — (2'X C — JxJy ~ JxJ Z )/2, 

r 4uy = (2Jg- y?- yf)/2, (8) 

T4UZ — {2J Z JzJx JzJy) /2j 

while the T 5u octupole operators are given by 

T 5ux = \/l5(Xjg - Xg)/6, 

T 5uy - Vl5(JyJl - JyJp/6, (9) 

T 5uz = VTE(J Z J% - J z J 2 )/6. 

For the evaluation of multipole susceptibility, we resort 
to the numerical renormalization group (NRG) method, 25 in 
which momentum space is logarithmically discretized to in- 
clude efficiently the conduction electrons near the Fermi en- 
ergy. In actual calculations, we introduce a cut-off A for the 
logarithmic discretization of the conduction band. Due to 
the limitation of computer resources, we keep m low-energy 
states. In this paper, we set A=5 and m=3000. Note that the 
temperature T is defined as T=K^ {N ^/ 2 in the NRG calcu- 
lation, where N is the number of the renormalization step. 

In Figs. l(a)-(e), we show the temperature dependence of 
T X T x for £ =5, £i=2, and E 2 =0.5 with A=1(T 5 . We note 
that such a small value of A is taken to reproduce the CEF 
level scheme for Pr-based filled skutterudites. 19 ' 20 Note also 
that we show one of Xx 7 ' s m eacn irreducible representation, 
since susceptibilities in the same irreducible representation 
take the same values owing to cubic symmetry. First of all, in 
the degenerate region with very small A, electron-hole sym- 
metry approximately holds within the j=5/2 sextet. Namely, 
the results for n=l and n=2 are quite similar to those for n=5 
and n=4, respectively, as observed in Figs. l(a)-(e). For n=\ 
and 5, at low temperatures, we observe that T4 U magnetic oc- 
tupole fluctuations are predominant. Note that the fluctuations 
of the T4 U dipole also survive, since it belongs to the same 
symmetry as the T4 U octupole. Also, for n=3, there remain 
magnetic fluctuations of the T4 U dipole and octupole. 

For n=2 and 4, on the other hand, we find significant fluc- 
tuations for the r 5g and r 3g quadrupoles, in addition to those 
of the T 4u dipole and octupole. Let us explain the results for 
the case of n=2. As shown in Fig. 1(f), for A>0, the local 
/-electron ground state is a T\ singlet composed of two r 7 
electrons. The first excited state is a T 4 triplet formed by T 7 
and Tg electrons, while the second excited state is a T 5 in- 
cluding a couple of T 8 electrons. Even if the Ti singlet is the 
ground state, there exists a T 4 triplet state with very small 
excitation energy. This T 4 triplet state carries a quadrupole 
moment. Thus, significant quadrupole fluctuations remain for 
n=2. The same explanation works also for n=4 due to the 
electron-hole symmetry for very small A. 
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Fig. I. Multipole susceptibilities vs. temperature at A=10 for (a) n=l, 
(b) n=2, (c) n=3, (d) n=4, and (e) n=5. (f) Local / 2 states in the j-j 
coupling scheme. 



Fig. 2. Multipole susceptibilities vs. temperature at A=0. 1 for (a) n=l, (b) 
n=2, (c) n=3, (d) n=4, and (e) n=5. (f) Local /-electron ground states for 
n=3, 4, and 5 with large A in the j-j coupling scheme. 



Now we consider the case with larger A, since the CEF po- 
tential is easily changed due to the substitution of rare-earth 
ion R, transition metal ion T, and pnictogen X for the filled 
skutterudite RT4X12. In Fij^s. 2 (a)-(e), we show the tem- 
perature dependence of T\x f° r A=0.1. For n=l, the low- 
temperature behavior of multipole susceptibilities is almost 
unchanged from those for A=10 -5 . Since one / electron is 
always located in a T7 orbital, Tg is irrelevant for the case 
of n=l. Since the Kondo temperature is very low for 1^=0.05, 
the local moment of the T7 electron still persists in the present 
temperature range, but it should be screened eventually for 
extremely low T. For n=2, the Ti singlet ground state is well 
separated from the T4 triplet excited state and thus, we do not 
observe any fluctuations at low temperatures. 

For ro=3 and 5, we note another electron-hole relation in 
the Tg quartet for very large A, as is easily understood from 
Fig. 2(f). In fact, the main results in Figs. 3(c) and (e) look 
similar to each other. In particular, we find that T2 U octupolar 
fluctuations are predominant at low temperatures for n=3 and 
5. Note that r.i u and T^ u octupolar fluctuations also become 
significant for n=5. Enhancement of octupolar fluctuations is 
a characteristic point for n=3 and 5 in the region of large A. 
Finally, for n=4 with large A, since the Tr levels are fully 
occupied, a couple of electrons in the Tg quartet form a local 
r$ triplet, leading to dominant magnetic fluctuations. 

In order to see the change of multipole fluctuations, it is 
convenient to show the Curie constant for multipole suscep- 
tibility, defined as C T ^=T X T x at T=l(r 10 . In Figs. 3(a)-(e). 
we plot C^ 7 vs. A for each n. In Fig. 3(f), we summarize 
the multipole properties. For n=l, as mentioned above, the 
multipole properties do not depend on A. Note again that 
for T<C10 -10 , the local moment of theT^ electron should be 
screened eventually. For n=2, in the region of A<5xl0~ 5 , 
both magnetic and quadrupole fluctuations are dominant, 



while for A>5xlCP 5 , multipole fluctuations are all "dead". 
For n=3, in the region of A<0.02, T4 U magnetic octupole 
and dipole fluctuations are dominant, even though the abso- 
lute values are small. On the other hand, for A>0.02, F2 U 
octupole fluctuations turn to be dominant. For n=4, T4 U dipo- 
lar fluctuations are always predominant, but quadrupolar fluc- 
tuations gradually die out as A increases. Finally, for n=5, 
in the region of A<2x 10~ 4 , T4 U magnetic octupole fluctua- 
tions are predominant, while for A>2x 10~ 4 , fluctuations of 
the T2u octupole are dominant and those of T4 U and Tsu are 
subordinate, similar to the results for n=3 with very large A, 
due to electron-hole symmetry in the r§ subspace. 

Next we discuss relevance of our results to RT4X 12 . Al- 
though the actual system is not exactly described by the im- 
purity model, the multipole properties are deduced from the 
dominant fluctuations. The case of n=\ would seem to cor- 
respond to the case of R=Ce, but in filled skutterudites, the 
cerium ion is considered to be tetravalent, rather than triva- 
lent. Thus, it may be difficult to claim a direct relation be- 
tween the n=l results and CeT4Xi2. For n=2, on the other 
hand, we expect a clear correspondence to the case of R=Pr, in 
which we confirm significant magnetic and quadrupole fluctu- 
ations at low temperatures. The range of A with magnetic and 
quadrupole fluctuations is very narrow, which may be related 
to the experimental fact that the appearance of exotic super- 
conductivity is very limited in Pr-based filled skutterudites. 

The cases of n=3 and 5 correspond to R=Nd and Sm, re- 
spectively. For both cases, in the region of small A, the pre- 
dominance of magnetic fluctuations suggests that magnetic 
ordering occurs in actual systems, although we have not de- 
termined whether it is ferro or antiferro. If it is possible 
to increase A experimentally for Nd- and Sm-based filled 
skutterudites, we expect exotic electronic properties to result 
from octupolar fluctuations, including T2 U and T5 U . In fact, 
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Fig. 3. Curie constant for each multipole susceptibility as a function of A 
for (a) n=l, (b) n=2, (c) n=3, (d) n=4, and (e) n=5. (f) Summary for the 
types of relevant multipole. Here "mag." denotes magnetic, while "quad." 
indicates quadrupole. The boundary between small and large A depends 
on n, but we simply ignore such a quantitative point in this table. 



future developments, (ii) We have suggested T2 U for the dom- 
inant type of octupole moment for n=5, but in the j-j coupling 
scheme, the effect of T n symmetry is not fully included. Since 
1?4 U and 1?5 U octupole moments should also occur in actual 
filled skutterudites, we cannot conclude which is the winner 
among octupole moments for R=Sm from the present results. 
To resolve this point, it is necessary to analyze a model in- 
cluding seven / orbitals, as has been done in Ref. 20. This is 
another future problem. 

In summary, we have analyzed the multiorbital Ander- 
son model to evaluate multipolar susceptibilities. It has been 
found that for n=l, magnetic fluctuations are always signifi- 
cant at low temperatures. For n=2, electronic states are dom- 
inated by magnetic and quadrupolar fluctuations at low tem- 
peratures, while multipolar fluctuations are dead for large A. 
For n=3 and 5, at small A magnetic fluctuations are dominant, 
while for large A, T2 U octupolar fluctuations are found to oc- 
cur. For n=4 at small A, the situation is quite similar to the 
case of n=2, but for increasing A, quadrupolar fluctuations 
gradually die out. 
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thanks R. E. Walstedt for critical reading of the manuscript. 
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of Science and by the Ministry of Education, Culture, Sports, 
Science, and Technology of Japan. The computation has been 
done using the facilities of the Supercomputer Center of In- 
stitute for Solid State Physics, University of Tokyo. 



for SmRu4Pi2, signs of octupolar ordering have been sug- 
gested by elastic constant measurements 26 and /iSR experi- 
ments. 27 When we turn our attention to actinides, quite re- 
cently, NpFe4Pi2 has been successfully synthesized. 28 Since 
the neptunium ion is considered to be tetravalent, Np-based 
filled skutterudites correspond to n=3. Thus, it may be natural 
to consider that the low-temperature properties of Np-based 
filled skutterudites are dominated by magnetic fluctuations. 

It is difficult to find actual filled skutterudite materials cor- 
responding to n=4, but Pu-based filled skutterudites may be 
good candidates if they can be synthesized, since actinide ions 
are tetravalent in the filled skutterudite structure. For n=4 with 
relatively small A, based on an approximate electron-hole re- 
lation in the j-j coupling scheme, one expects to find elec- 
tronic properties similar to those of Pr-based filled skutteru- 
dites. In particular, exotic superconductivity may also occur 
in Pu-based filled skutterudites with relatively high transition 
temperatures, since hybridization is large in actinide com- 
pounds owing to the difference in the nature of 4/ and 5/ 
electrons. This point should be carefully discussed, when Pu- 
based filled skutterudites are synthesized at some future point. 

Two comments are in order, (i) In this paper, we have sim- 
ply ignored the effect of the rattling motion of rare-earth ion, 7 
but it is interesting to see whether magnetic and quadrupolar 
fluctuations are enhanced or not due to the effect of rattling 
in Pr-based filled skutterudites. This point can be investigated 
based on an extended model incorporating further electron- 
phonon interaction terms into Hi oc . This is among possible 
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